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Abstract: The residual dependence index of bivariate Gaussian distributions is determined by the 
correlation coefficient. This tail index is of certain statistical importance when extremes and related rare 
'. events of bivariate samples with asymptotic independent components are being modeled. In this paper 

^jr^' we calculate the partial residual dependence indices of a multivariate elliptical random vector assuming 

^ , that the associated random radius is in the Gumbel max-domain of attraction. Furthermore, we discuss 

' the estimation of these indices when the associated random radius possesses a WeibuU-tail distribution. 

■ Key words and phrases: Partial residual dependence index; Gumbel max-domain of attraction; Weibull- 

tail distribution; Elliptical distribution; Quadratic programming problem. 



1 Introduction 



»^ ' Let {Xi,X2) be a bivariate elliptical random vector with stochastic representation 



(Xi,X2) = R(Ui,pUi + ^1- p^U2), pe{-l,l), (1.1) 



where the positive random radius R is independent of {Ui, U2) which is imiformly distributed on the unit 

^ ■ circle oilR^. Here = stands for equality of distribution functions. A canonical example of a bivariate 

i elliptical random vector is when is Chi-square distributed, which implies that Xi , X2 are standard 

10 i Gaussian random variables with mean 0, variance 1 and correlation coefficient p := E{XiX2}- It is 

well-known (see e.g., Reiss and Thomas (2007)) that in the Gaussian model the correlation coefficient p 
does not influence the asymptotic dependence of the components. Roughly speaking this means that the 
sample extremes of Gaussian random vectors are asymptotically independent. 

A tractable extension of the Gaussian model is the elliptical one, where R is some general positive random 



I variable with distribution function F. By Lemma 12.1.2 of Berman (1992) 



Xi=X2 = RUi = RU2 (1.2) 



implying that distribution function of Xi (denoted below by Q) is continuous. Clearly, the joint depen- 
\ dence function of Xi and X2 is influenced by p. 

Several authors have considered elliptical distributions for modelling of rare events with specific appli- 
cations in insurance and finance. Recent contributions in this directions are Peng (2008), Li and Peng 
(2009). In the first mentioned article the author deals with the situation that the marginal distributions 
are regularly varying implying that the components are asymptotically dependence. In such a model, 
also considered in Kliippelberg et al. (2007), novel estimation techniques are presented. 
In Berman (1983,1992), Hashorva (2005a,b), Abdous et al. (2005,2008) further probabilistic results are 
obtained when F is in the Gumbel max-domain of attraction which implies that the components Xi and 
X2 are asymptotically independent. 

In this paper, with motivation from above mentioned contributions, we focus on the quantification of the 
asymptotic dependence of elliptical random vectors. 

An interesting measure of the asymptotics dependence (see e.g., Peng (1998, 2007), de Haan and Ferreira 
(2006), Reiss and Thomas (2007)) is the function x(m) defined by 

P{Xi >u,X2>u} ^ 
^(")^= P{X^>u} 



2 



If for some constant c e [0, 1] we have 

lim x{u) = ce (0,1], (1.3) 

u — ►oo 

then Xi and X2 are said to be asymptotically dependent. In our setup of bivariate elliptical random 
vectors with stochastic representation (jl.ip this is the case when R has distribution function F in the 
Frechet max-domain of attraction (or equivalently, F is regularly varying with positive index 7). See 
Berman (1992) or Hashorva (2005a, 2006b) for further details. Important statistical applications can be 
found in Kliippelberg et al. (2007). 

In both other cases of max-domain of attraction, i.e., F is in the Gumbel or the WeibuU max-domain 
of attraction we have (see Hashorva (2005a)) c = 0, which means that Xi and X2 are asymptotically 
independent. 

In extreme value theory asymptotic independence is a nice property, however, c = in ()1.3I) merely means 
that P{Xi > u, X2 > u} converges faster to than the marginal survival probability P{Xi > u} (if 
u — > 00). 

One successful approach to model the asymptotic independence is the estimation of the residual depen- 
dence index 77 e (0,1] (see e.g., Peng (1998, 2007, 2008), de Haan and Peng (1998), or de Haan and 
Ferreira (2006)). We note in passing that recent ideas in testing asymptotic can be found in Hiisler and 
Li (2009). Now, information about rj is available if for any x,y positive 

S^ix,y):=^^^ ^ 5(x-,2/)e (0,00), 00, (1.4) 

'^'11(15 1 ) 

with 

Su{x, y) := P{Q{X,) > 1 - x/u, Q{X2) > 1 - y/u}, u>0, 
since for any c > and for some 77 e (0, 1] we have the important scaling relation 

S{cx,cy) = c^/'^S{x,y). 

Furthermore, the fimction ^^(l, 1) is regularly varying at infinity with index —1/rj. Other authors refer 
to ri as the coefficient of tail dependence (sec e.g., Resnick (2002), or Reiss and Thomas (2007)). 
In this paper we consider the problem of calculating the residual dependence index r/ for the bivariate 
random vector [Xi,X2) with stochastic representation (jl.ip assuming that the distribution function F 
is in the Gumbel max-domain of attraction. We show that r] does not always exist. In certain instances 
when it exists we prove that -q is defined in terms of p and the WeibuU tail-coefficient 9 (see below (j2.11|) ). 
In Section 3 we propose an estimator of the residual dependence index rj. Definition, calculation and 
estimation of the partial residual dependence index for multivariate elliptical distributions are placed in 
Section 4. In the multivariate setup the partial residual dependence indexes (if they exit) are determined 
by the unique solution of specific quadratic programming problem, and the WeibuU tail-coefficient 9. 
Proofs of all the results are relegated to Section 5 (last one) . 



2 Calculation of the Residual Dependence Index 

Let (Xi, X2) be an elliptical random vector with stochastic representation (|l.ip . and let R be the positive 
associated random radius with distribution function F. We assume in the following F(0) = and 
F{x) < l,Va; > 0. If Xi and X2 are standard Gaussian random variables, then it is well-known that 
(see e.g., Reiss and Thomas (2007)) Xi and X2 arc asymptotically independent for any p € (—1,1). 
Furthermore, the residual dependence index 77 is given by 

r;:= (1 + p)/2 G (0, 1) 

and (sec p. 322 in Reiss and Thomas (2007)) 

5„(1,1) = (1 + 0(1)) ^"f/" {'^n)-p'^^+p\\nu)-p'^^+P\-^'(^+P\ u^oo. (2.5) 
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In our notation o(l) means liniu_,oo o(l) = 0. 

When Xi and X2 are independent, then i?^ = + X| is Chi-squared distributed with 2 degrees of 
freedom. The distribution function of i? is in this case in the max-domain of attraction of the Gumbel 
distribution A(a;) = exp(— exp(— a;)), a; G 5?. From the extreme value theory we know (see e.g., Resnick 
(1987), Reiss (1989), Falk et ah (2004), Embrechts et al. (1997), or de Haan and Ferreira (2006)) that the 
distribution function F of i? is in the max-domain of attraction of A, if for some positive scahng function 
w we have 

i — F(u + x/w(u)) / X w ^ „N 

Hm ^ — ^ " = exp(-a;), eM. (2.6) 

ti— >oo 1 — F[u) 

In the Gaussian case 77 is strictly less than 1 and w{u) = (1 + o{l))u, u > 0. In the more general elliptical 
setup of this paper it turns out that interesting cases for calculation of are when 

w{u) = u''-^L{u), ee[0,oo), (2.7) 

with L a positive slowly varying function at infinity satisfying lim„_too L{cu) / L{u) = l,Vc > 0. We refer 
to 6 as the Weibull tail-coefficient index (see Girard (2004)). 

We show below that the elliptical model exhibits two main two main features, namely: a) the residual 
dependence index rj (when it exists) depends on both p and 9, being in fact an increasing function of p 
and \/9, and b) it is possibly to have rj = 1 when 6 = Q and lim„_+oo L{u) = 00. More interestingly, both 
X, Y are asymptotically independent even when ry = 1 . 
The main result of this section is the following theorem. 



Theorem 2.1. Let {Xi^ X2), p G (—1, 1) be a bivariate elliptical random vector with stochastic represen- 
tation (fTT]) . A ssume that R has distribution function F which satisfies p.6p with some positive scaling 
function w. 
(i) Suppose that 

lim = 00, with ap := ^2/{l + p) > 1 (2.8) 

u^oo W[U) 

holds, then 

lim Su{x,y) = 00, ifx>l,y>l, lim S'„(a;, ?/) = 0, ifx,ye{0,l). (2.9) 
(a) If for some 9 £ [0, C5o) 

lim = VOO, (2.10) 

u—foo w{u) 

then for any x,y ^ (0, oo) 

61/2 



hm Su{x,y) = (xy)i/(2.), -.J^/] 



e(0,i], (2.11) 



and 5*^(1, 1) is regularly varying at infinity with index —l/rj. 

(Hi) Let denote the inverse of the distribution function of Xi. As u ~* oo we have the asymptotic 
expansion 

5„(1,1) = (l + o(l)) ^ . )*)'' b^{u):^a,Q'\l^l/u). (2.12) 

27r(l — p)^ b^,(u)w[b^,(u)) 

Remarks 2.2. 1) The scaling function w in (j2.6p can be defined asymptotically by (see e.g., Resnick 
(1987)) 

, , (1 + o(l))[l - F{u)\ 

^(.) := ' ^^\,\^;,^ — oo. (2.13) 
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Further, we have 

lim uw{u) = oo. (2-14) 

u — ^oo 

Hence in the model (j2.7p the Weihull tail- coefficient is necessarily non-negative, and if 9 ~ 0, then we 
need to suppose further that limti_>oo L{u) = 0. Two interesting distributions with 9 = and L{u) = 
clnUjC e (0,oo),u > are Benktakder type I and Lognormal one, see Embrechts et al. (1997), pp. 149- 
150. 

The interesting model 9 = is kindly suggested by the referee of the paper. A striking feature of this 
model is that residual dependence index rj equals 1, thus not depending on p at all. Furthermore, in view 
of Hashorva (2005a) X and Y are asymptotically independent. 

2) In view of (|1.2p Xi is a product of R and an independent random variable Ui. Asymptotics of random 
products are investigated by several authors, see for recent results Tang and Tsitsiashvili (2003,2004), 
Tang (2006,2008). 

Next, we present three examples. 

Example 1. Let {Xi^X^), p be as in with associated random radius i? ~ A. Clearly, the unit 

Gumbcl distribution A is in the Gumbel max-domain of attraction. An admissible choice for the scaling 
function is w{u) = l,Vw > 0. Consequently, (|2.10p holds with = 1, implying that 1) is regularly 

varying with index —(0.5 + p/2)~^/^. 



Example 2. Under the setup of Example 1 we assume further that R has distribution function F in the 
Gumbel max-domain of attraction with the scaling function w{u) = exp(au), m > 0, with a some positive 
constant. Such F exists and can easily be constructed if we assume that F possesses a density function 
/, requiring further /(it)/[l — -F(u)] = w{u),\/u > 0. For this choice of the scaling function w (|2.8p holds. 
Hence lim„^oo Su^x, y) = for any x,y <E (0, 1). 



Example 3. [Kotz Type III] Again with the setup of Example 1 if for all large u 

P{R>u} = (1 + o(l))A'w^exp(-™''), K > 0,9 > 0,N eM, 



(2.15) 



then we refer to {Xi, X2) as a Kotz Type III elliptical random vector. Since we assume that 9 is positive, 
R has distribution function F in the Gumbel max-domain of attraction with the scaling function 

w{u) = {1 + o{l))r9u'^^'^ , u^oo. 

Consequently, (|2.10p holds and ry = £ (0, 1]. Next, if we define b{u) := Q~^{1 - l/u),u > with Q 
the distribution function of Xi , then Theorem 12.11 implies 



Su{l,l) 



In view of Theorem 12.3.1 in Berman (1992) and the fact that Q is symmetric about (recall (1.1)) 



l-Q{u) = -P{Xl>u} 



-P{R^Ui >u} = (1 + 0(1)) 



K 



V2^ 



i^~''/2exp(-ra''), u-^oo. 



where Ui is beta distributed with parameters 1/2,1/2 being independent of R. Hence we may define 
h{u) asymptotically as (see Embrechts et al. (1997)) 



b{u) = [r^^lnuf 



1 



(1 + 0(1)) 



Olnu 



[N - 9/2) ln(r-i In u) /O + lnK-]^ ln(27rr6 



Thus we arrive at (set A := a^) 

(l+o(l))^^ ' 



Su{lA) 



K{l-pf 



\2t:9 J 



(Inw) 



{l~X)N/e+\/2~l^-X 



u — > 00. 
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In the special case 

K^l, r = l/2, 9 = 2, N = Q, ap = ^2/{l + p), 
which holds in particular if both Xi and X2 are standard Gaussian random variables we retrieve (|2.5p . 

3 Estimation of 77 in the Weibull Model 

In view of Theorem 12. II if the scaling function w is regularly varying with index 9—1, then the residual 
dependence index rj is defined in terms of p and 9. Let {Xki, Xk2), k = 1, . . . , n be a sample of bivariate 
elliptical random vectors with stochastic representation (|l.ip (where p € (—1,1) is assumed). Then a 
non-parametric estimator p„ of p is given by (see e.g. Peng (2008), Li and Peng (2009)) 

Pn sin(7rf„/2), n> 1, (3.16) 

where f„ is the empirical estimator of the Kendall's tau. 

Good performing estimators of the so-called Weibull tail-coefficient are the Girard and Zipf estimators, 
see e.g., Girard (2004). Referring to the aforementioned paper we say that the random radius R ^ F 
possesses a Weibull-tail distribution if 

1 - F{x) = cxp{-H{x)), H-^ix) = M{t : H{t) > x} = x^/^Li^x) (3.17) 

holds with Li a positive slowly varying function at infinity. Gardes and Girard (2006) and Diebolt et al. 
(2008) give several examples of Weibull-tail distributions. Prominent instances are the Gaussian, Gamma 
and extended Weibull distributions. 

By the properties of slowly varying functions (see e.g., de Haan and Ferreira (2006)) we may write (|3.17p 
alternatively as 

1-F{x) = exp(-.T^L2(a;)), (3.18) 
where L2 is another slowly varying function which is asymptotically unique. 

From the estimation point of view, a tractable class of the Weibull-tail distributions is constructed when 
F is in the Gumbel max-domain of attraction with the scaling function w defined asymptotically by 

wiu) = TTTTTT' ce(0,oo), (3.19) 
1 + ti{u) 

where ti is a regularly varying function at infinity with index k := 9p, p £ (—00, 0) implying 

1-F{u) = cxp{-cu\l+t2{u))), u>0, (3.20) 

where <2 is another regularly varying function at infinity with index n. 

Under the assumption (|3.19p it follows that (see Berman (1992) or Hashorva (2005a)) 

P{Xi > u} = cxp{-cu\l+t3{u))), ueM, 

with ^3 again a regularly varying function at infinity with index k. 

Next, assume that the associated random radius i? defining the random sample {Xki, Xk2), k = 1, . . . ,n,n > 
1 possesses a Weibull-tail distribution F such that (|3.19p holds. Write Xi-n < ■ ■ ■ < Xn-.n for the associ- 
ated order statistics of Xn, . . . ,Xni- Following Gardes and Girard (2006) we might estimate 9 by 

11" 

2 — 1 

with 1 < A^yj < n, > 0, n > 1 given constants satisfying 

lim kn — C50, lim — — 0, lim log(T„/A:„) ~ 1, lim \/hnb{[og{n/kn)) A G-/R, 

n—*oo n— ^ctt Ti n — >-oo n— i-oo 
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where the function 6 is a regularly varying function with index 9 appearing in a second order asymptotic 
condition imposed on F (being thus related to Li). 

Asymptotic properties of On are discussed in Gardes and Girard (2006) and Dicbolt et al. (2008). Next, 
based on our main result we propose an estimator for the residual dependence index i] given by 

fin ((l + p„)/2)'"^', n>l. (3.21) 

Asymptotic properties of ^„ follow by utilising the asymptotic properties of both /5„ and 0„ . 
We note in passing that the constant c can be estimated by 

1 ^ log(n/») ^ 

•^n T n — 'i+ 1 : n 

4 Partial Residual Dependence Index 

Consider X := {Xi, . . . ,Xfc)^, fc > 2 an elliptical random vector mIR'' with stochastic representation 

X = RA^U, (4.23) 

where R is again the positive associated random radius of X with distribution function F independent 
of t/ := {Ui, . . . ,Uk)^ which is uniformly distributed on the unit sphere oiM^ and A ^W&^^ is a non- 
singular matrix (here ^ stands for the transpose sign). The distribution function of the random vector X 
is determined by the positive definite matrix S := A^ A, the distribution function F and the distribution 
function of U (which is known). See Cambanis et al. (1981) or Fang et al. (1990) for more details on 
elliptical distributions. 

We assume in the following that F(0) = 0, F(a;) < l,Va; > and S is a correlation matrix i.e., all the 
entries of the main diagonal equal 1. If the distribution function F is in the Gumbel max-domain of 
attraction, then each pair Xi^Xj^i ^ j < k is asymptotically independent. If further the scaling 
function w satisfies (|2.10p . then by Lemma 12.1.2 in Berman (1992), Proposition 3.4 in Hashorva (2005a) 
and Theorem 12. H it follows that the residual dependence index rjij (related to {Xi,Xj)) is 



with pij G (—1,1) the ij-th entry of E. 

Let / be a given non-empty index subset of {1, . 

2, . . . , fc. Next, define Su.i{x) by 



. . , fc}. By the assumption on S we have Xi = Xi, i = 



Suj{x) := p{q{X,) > 1 - ^,V^ el}, u>0,x {xi, . . .,XkV e (0,c»)^ 
with Q the distribution function of Xi. 

If 'S'ti,7(l) is regularly varying with index —1/77/,?// e (0, 1]. then we refer to 777 as the partial residual 
dependence index of the subvector Xj := {Xi,i G /)^, or shortly as the partial residual dependence 
index. 

The submatrix of E obtained by retaining the rows and the columns of S with indices in J and /, 
respectively, (assume / has less than fc elements) is denoted by 'Sji,J := {l,...,fc} \ /. We define 
similarly E77 and xj for x elR^. 

Since in our model S = A'^ A is positive definite the inverse matrix of S77 exits. Next, we write ai 
for the unique solution of the quadratic programming problem 

minimise the objective function yjY.Jjyj, y := (j/i, . . . , j/^.)^ G 5?'^, j/i > 1, yi<El. (4.24) 

In the multivariate setup again two interesting features are observed, namely a) the partial residual index 
rji depends on both ai and 6 (when it exists), and b) when 9 = 1, then rji = 1 for all / C {1, . . . , fc}. 
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Theorem 4.1. Let X be an elliptical random vector inW&^k > 2, with stochastic representation (|4.23p . 
Assume that the associated random radius R is almost surely positive with distribution function F satis- 
fying (j2.6p . // the scaling function w satisfies (j2.10p . then for any non-empty index set I C {1, . . . ,k} 
with m < k elements and any x G (0, oo)'^ we have 

\un^^ = (n^?)' 7, :=«rHeJS^i^li^)e(0,oo), e K, (4.25) 

'^uJ\^) jeK 

where K is a unique subset of I with I > elements such that Y,J^\^1k is a vector with positive ele- 
ments and ej is the j-th unit vector inlRl'. Furthermore, if M :^ I \ K is not empty, then the vector 
I]ifM5^j\fj\,f Ia/^Ik has non-negative components and S'„_/(l) is regularly varying with index —l/rjj where 

m ■■= «7''e ((i7s]7/i,)"'',i). (4.26) 

Remarks 4.2. 1) Estimation of the partial residual dependence index rji related to a given index set I 
requires estimation of the attained minimum aj of the related quadratic programming problem and the 
Weibull tail- coefficient 6. For any i,j,i ^ j an estimator of pij (the ij-th entry of Y,) can be defined 
analogously to p.l6p by 

Pi3,n ■= sin(7r'ry_„/2), n > 1, 
with fij^n the corresponding empirical estimator of the Kendall 's tau. 

An estimator of a/ can be constructed if we already have estimated the precision matrix E^/ • Estimation 
of TjJ^ is recently discussed for Kotz distributions in Sarr and Gupta (2008). 

If Oil, n denotes an estimator of aj, and 0„ an estimator of the Weibull tail-coefficient, then in view of 
our results we can estimate the partial residual index r]j by 

m,n a7>- (4-27) 

2) In the case that (|2.8p holds with aj instead of ctp, then we cannot define rjj. 

3) IfEj^lj has positive elements, then the subset K in Theorem \4.1\ equals I. This is in particular the 
case if I has only two elements, or when the non-diagonal elements ofT. are all equal, say to p G (^li 1)- 
If K I, then for estimating aj, we also need to identify the elements of K , which is not an easy task 
in general. 

4) It is well-known that the solution of the attained minimum aj of the quadratic programming problem 
above is related to the exact tail asymptotics of the Gaussian random vectors, see Dai and Mukherjea 
(2001), or Hashorva (2005b, 2007b) for more details. 

5) In the particular case w is given by (j2.7p with 6 = 0, then 77/ = 1 for all subsets I of {1, . . . , k}. 
Furthermore, E does not influence rjj, which is in particular the case for multivariate lognormal distri- 
butions. 

Next we consider the trivariate setup in some details. The foUowing lemma gives an explicit formula for 
a I, I = {1,2,3}, which is useful for the estimation of aj. 

Lemma 4.3. Let E G IR,^^^ be a positive definite correlation matrix (with I's in the main diagonal) 
and non-diagonal entries pij G ( — l,l),i 7^ < 3. Define pmin iiiin(pi2, P13, /923) md set a := 

niina;,>i,i=i,2,3 x^E~'^x. 

(i) // 1 + 2p„ii„ - pi2 - Pis - P23 > 0, then we have (here 1 = (1, 1, 1)^) 

a = 

^ 3 - 2(pi2 + pi3 + P23) - P12 - Pi3 - P23 + 2(pi2Pi3 + P12P23 + P13P23) 28) 

1 + 2pi2Pi3P23 - P?2 - Pl3 - P23 

(ii) If 1 + 'i.pmin — P12 — Pi3 ~ P23 < 0, then there exists a unique index set C {1, 2, 3} such that 

Pmin = Pij < mm pik. (4.29) 
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Moreover we have 

a = llE],'^lK^{l,lV^Aihl) = -^. (4.30) 

Example 4. [Kotz Type III, 3-dimensional Case]. Let X be an elliptical random vector in J?'^ with 
stochastic representation (|4.23p . where the matrix A is non-singular and set E := A. We denote by 
Pij the ij-th entry of S. Assume that pa = l,i = 1, . . . , k and R satisfies (|2.15p as u ^ oo. Again we 
refer to X as a Kotz Type III random vector. In view of Lemma 12.1.2 in Bcrman (1992) we have for 
any index set / = {k, 1} C {1, 23} with two elements 

iXk,Xi) = R(Ui,pkiUi + ^\-pIiU2), 

where {Ui,U2) with uniform distribution on the unit circle of IR^ is independent of R. Hence we can 
estimate pki as in p.l6p . Let Pi2,m Pis.m P23,m'n > 1 denote these estimators. 

Next, consider the case / = {1,2,3}. In view of Theorem 14.11 and Lemma 77 r = with a defined 
in gIMl). If 

1 + 2 min(pi2,„, PlZ,n, P23,n) ~ Pl2,n " Pl3,n ^ P22,n > 0, 

then the estimator of a is obtained by plugging in the estimators pi2,n, Pi3,n, P23,n- Otherwise, we estimate 

Pmin.n ■= niin(/5i2.„ , /5i3,„ , P23,n) , n > 1, 

and obtain the estimator of a be plugging in Pmin.n in (|4.30p . The WcibuU tail-coefficient 9 can then be 
further estimated as previously discussed in Section 3. 



5 Proofs 

Proof of Theorem 12.11 Let Q be the distribution function of Xi with inverse {Q is a continuous 
function, see e.g., Berman (1992)). The Gumbcl max-domain of attraction assumption on F implies (see 
e.g., Reiss (1989), or de Haan and Ferreira (2006)) 

w{b{u))[Q-^{l-x/u)-b{u)] -Inx, u ^ oo (5.31) 

locally uniformly for x £ (0,oo), with b{u) := Q^^{1 — l/u),u > 0. Next set 



w*{u) ■.= w{apb{u)),u>0, Up := yj2/{l + p) > 1, pe(-l,l). 
For any u, x, y positive we may further write (recall Xi = X2) 
Su{x,y) = p{q{X^)>1--M{X2)>1--] 

L U U J 

= p|Xi >Q-i(l--),X2 >Q-i(l--)| 
L u u J 

= P{X, > bin) - (1 + oil))^^ X2 > bin) - (1 -f 0(1))^^}. 
In view of Theorem 5 in Hashorva (2007a) for any s, t positive 

lim pIw*{u)(Xi - b(u)) > s,w*(u)(X2 - b(u)) > t Xi > b(u),X2 > b(u)] = P{X[ > s,X!,> t}, 



holds with X[ , X2 two independent exponentially distributed random variables with mean Ap := ■\/2(l + p). 
Hence \i x,y £ (0, 1), then — ln.x, — In y e (0, 00), thus (|2.9p follows easily. For any a; > 1 and ?/ > 1 we 
may write 

lim Su{x,y) = lim — — — — — = 00. 
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Next, if ((TTUl) holds, then 

w*(u) 



G-\l-x/u)~b{u) 



I, u ^ oo 



locally uniformly for any x > 0. Consequently, with the same arguments as above for any x,y Cz (0, 1] we 
obtain 

lim Suix,y) = P{X[ > -al'^ \nx, X!^ > -a^-Hny} 

= exp(^^^ln(xy)j = exp(^^ ln(a-2/) j =: S{x,y). 



The result for x S (1, oo) and y positive, or x positive and y G (1, oo) as well as the statement (iii) can 
be now established using directly Theorem 2 in the aforementioned paper. Since for any c, x, y positive 

S{cx,cy) = 5(a;,?/)c^/'', 

with 

/, ^ 61/2 

'?:=<=(^) e(0,l] 

the result follows. □ 
Proof of Theorem 14.11 Let / be a no-empty subset of {1, ... , k} with m < k elements. The random 
vector X] := {Xi,i G /)^ is again an elliptical random vector with stochastic representation (Cambanis 
et al. (1981)) 

Xi = RiBV, 

with positive associate random radius Rj, square matrix B such that B = E// and Y uniformly 
distributed on the unit sphere ofiR"' being independent of Rj. As shown in Berman (1992) the associated 
random radius Rj has distribution function Fj in the Gumbel max-domain of attraction with the same 
scaling fimction w as the distribution function of R. By Proposition 2.1 in Hashorva (2005b) there 
exist a unique subset K d I with I > elements such that 

°" ^w,m 1 y^^lh = 'LI^kk'^k > 0, 

y(^lH.^^ ,yi > 1,2— 1 , . . .,?n 

E^^lx has non-negative components, and if M := I \ K is not empty, then Y.km'^^im'^m ~ has 
non-negative elements. 

As in the proof of Theorem 12. II for any x G (0,00)*^ applying further Theorem 3.4 in Hashorva (2007b) 
we obtain 

Sujjx) _ ( Y^^l^ja",-^ 

with := ejT.j}j^lK > 0, and ej the j-th unit vector inJR'. We have further 



li^i^ = iU-r' =:^(^)G(0,oo), 



jeK j&K 

hence for any c > and any x = (xi, . . . , Xfc)^ G (0, oo)'^ we may write 

Sicxi,...,cxk) ^ (TT'^^'"'') 
S{xi, . . . , Xk) jeK 

Consequently, 777 = aj^ , thus the result follows. □ 
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Proof of Lemma [4.31 The proof of the first statement is shown in Lemma 3.2 in Hashorva and Hiisler 
(2002). We show next the second statement. Assume therefore that 

1 + 2p,„j„ - pi2 - Pl3 - P23 < 0, and p,nin = Pl2- 

Since 1 — max(p23,Pi3) > 0, then pi2 ~ P23 or pi2 — pis is not possible, hence pmin = Pi2- In view of 
the aforementioned lemma a = 2/(1 + P12); thus the result follows. □ 
Acknowledgement: I am in debt to the referee for deep professional review, very kind report, several 
suggestions and corrections. 
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